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Abstract. We have analysed the Rhodes/HartRAO survey at 2326 MHz and derived the global angular power spectrum of 
Galactic continuum emission. In order to measure the angular power spectrum of the diffuse component, point sources were 
removed from the map by median filtering. A least-square fit to the angular power spectrum of the entire survey with a power 
law spectrum Ci oc l~ a , gives a — 2.43 ± 0.01 in the I range 2 — 100. The angular power spectrum of radio emission appears 
to steepen at high Galactic latitudes and for observed regions with |6| > 20°, the fitted spectral index is a = 2.92 ± 0.07. We 
have extrapolated this result to 30 GHz (the lowest frequency channel of Planck) and estimate that no significant contribution 
to the sky temperature fluctuation is likely to come from synchrotron at degree-angular scales. 
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1. Introduction 

The accuracy with which we are able to derive the cosmolog- 
ical parameters from the angular power spectrum of the cos- 
mic microwave background (CMB) depends on the properties 
of Galactic and extra-galactic foreground emission. The future 
satellite missions Planck (ESA) and MAP (NASA) designed to 
perform high precision observations of the CMB have a wide 
frequency coverage in order to disentangle the various com- 
ponents of the microwave sky and provide clean CMB maps. 
The two satellites, therefore, will also provide new high qual- 
i ty data on Galactic dif fuse emission and extra-galactic sources 
(Bersanelli et al. 1996). However before these observations be- 



come available it is important to improve our knowledge of 
the statistical properties of the foreground emission in order to 
quantify the level of contamination of the CMB signal expected 
in the various channels and at the different angular scales. This 
will allow more realistic simulations of the missions to be per- 
formed and the analysis techniques to be refined. 

Below 50 GHz measurements of the CMB anisotropy are 
mostly affected by Galactic radio emission: a combination of 
thermal bremsstrahlung from ionized gas, synchrotron radia- 
tion from cosmic ra y electrons and possibly em ission from 
spinnin g dust grains (|de Oliveira-Costa et al. 1999 ). As pointed 
out by Tegmark et al. (2000) when trying to recover the an- 
gular power spectrum (APS) of the CMB, the uncertainties in 
the shape of the APS of a foreground component is a greater 
source of error than the uncertainties in the frequency depen- 



dence. The only available estimates of the APS of synchrotron 
emission comes from the analysis of the 408-MHz survey of 



Haslam et al. (1982) and 1420-MHz survey of Reich & Reich 



(1986). They suggest a power-law spectrum with a spectral in- 
dex around 2.5—3.0 down to the maps' resolution limit ~ 1° 



( [Tegmark & Efstathiou 1996| ; [Bouchet et al. 1996| ; Bouchet & 
Gispert 1999). 

Recently the Rhodes/HartRAO radio continuum survey at 
2326 MHz (Jonas et al. 1998) has become available. In this 
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paper we use these new data to derive the APS of the radio 
continuum at 2.3 GHz for all the observed sky and for high 
Galactic latitude regions. We then extrapolate the high Galactic 
latitude results to predict the level of contamination that might 
be expected from diffuse synchrotron emission in the lowest 
frequency channel of Planck. 

The paper is organized as follows. In Sect. ^ the data from 
the Rhodes/HartRAO survey are presented and the APS is de- 
rived for different regions of the sky, after point source sub- 
traction. In Sect. H the simulations which were performed to 
assess the reliability of our results are described. Simulations 
and results are discussed in Sect. || 

2. The angular power spectrum of radio 
emission at 2.3 GHz 

The Rhodes/HartRAO radio continuum survey at 2326 MHz is 
a single-dish s urvey of the sky w ith a resolution (at FWHM) 
of 20 arcmin ( Jonas et al. 1998 ). The survey made from the 
Hartebeesthoek Radio Astronomy Observatory (HartRAO) in 
South Africa covers the declination range —83° < S < 32° 
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for RA 360° - 180° and -83° < 8 < 13° for most of RA 
180° - 0° (in the RA range 90° - 15°, -80° < S < 13°). 
The observations are spatially oversampled both in RA direc- 
tion (6 arcmin) and in declination direction (1.5 arcmin). The 
uncertainty in the temperature scale is less than 5 per cent and 
the error in the absolute zero level is less than 80 mK in any 
direction. 

In order to analyse the global properties of the radio 
continuum emission at 2326 MHz, the raw data of the 
Rhodes/HartRAO survey were re-gridded into a HEALPix tes- 
sellation. HEALPix is a Hierarchical, Equal Area and iso- 
Latitude Pixelisation of the sphere designed to support effi- 
ciently: local operations on the pixel set, a hierarchical tree 
structure for multi-resolution ap plications and the g lobal Fast 
Spherical Harmonic transform (Gorski et al. 1999). This last 
property is crucial for fast computation of the APS of a field 
defined on the sphere. The recent CMB maps obtained from 
the observations of the balloon-born experiment Boomerang 
have also been pixelised according to the HEALPix scheme 
dde Bernardis et al. 20001). 



In the HEALPix tessellation the angular size of the pixels 
is determined by the nside parameter. For the Rhodes/HartRAO 
survey we used a HEALPix tessellation with nside=5\2 which 
corresponds to pixels of linear size 6.9 arcmin and there- 
fore provides an adequate re-sampling of the survey. The 
Mollweide projection of the raw data of the survey re-gridded 
into HEALPix is displayed in Fig. |l|. A great number of point 
sources are clearly visible in the map and the survey appears to 
be heavily affected by noise at the pixel scale. 

The APS of a field on a rectangular-domain is given by 
the Fourier components of the 2-point correlation function. The 
analogue in spherical geometry are the C/ components of the 
Legendre polynomial expansion 

C{9) = j-53(2I + l)C,fi(cos0) 
i 

If T(8, 4>) = 2~2t m a imYi m (9, 4>) is the spherical harmonic de- 
composition of the field defined on the sphere then the C\ co- 
efficients are simply given by: 
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The APS of the Rhodes/HartRAO survey has been ob- 
tained by performing a spherical harmonic analysis of the field 
in the HEALPix grid with the ANAFAST program (the soft- 
ware is part of the HEALPix package available at the web site 
http : / / www . eso . org/kgorski/healpix/). The re- 
sult of the analysis for the whole set of raw data are shown in 
Fig. ||, as well as the APS derived after imposing a Galactic 
cut with |b cut | = 10° and the APS derived for |b cut | = 20° []. 
The presence of an equatorial cut across the sphere affects the 
harmonic decomposition of the field. Performing a cut in the 
space domain is equivalent to multiplying the observed field by 



1 A Galactic cut with b — b cut , means that all the data in the latitude 
range — b cut < b < b cut , have been set to zero 



a window function and, in the frequency domain, this is equiv- 
alent to the convolution of the Spherical Transform of the field 
with the transform of the window function. As a Galactic cut is 
a discontinuous window function, its Spherical Transform is af- 
fected by aliasing. The resulting APS is affected by slow oscil- 
lations, which reflect the envelope of the transformed window 
function (which in this case is the combination of the incom- 
plete sky coverage with the Galactic cut), and fast oscillations 
due to the aliasing signal. The oscillations are clearly visible in 
Fig.| 

The cylindrical symmetry of a Galactic cut manifests itself 
as a "coherent" set of terms in m=0, therefore the effect of this 
window function on the measured APS can be mitigated by 
implementing a new definition of the APS in which one sets 
the a;o terms to zero, that is: 
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Eq. |2]is the definition (which we shall refer to as Modified APS, 
hereafter) that has been used to obtain the power spectra from 
the data with a Galactic cut in Fig. ^]and Fig. ^ (grey-lines). 

The window function of a Gaussian beam with dispersion 
Oh is given by: 

Wi«exp[-iofl(i + l)] 

which means that for a beam with a FWHM of 20 arcmin 
the beam suppression of the sky APS starts becoming ap- 
preciable (transmission factor < 0.9) at I ~ 200 (<r& = 
FWHM/V81n2). We model the APS of this component of the 
radio sky as a single power law: 

Ci = A 2 r a 

Up to I — 200, a linear least-square fit to the APS of all raw 
data (black line in Fig. §) gives A = 0.63 ± 0.14 K and a = 
2.11±0.06. The APS at high Galactic latitude cannot clearly be 
described by a single power law up to I — 200, however if one 
attempts a linear fit up to I = 100, where the curves flatten, one 
obtains: A = 0.2i0.1K and a = 2.38 ± 0.05, for |b| > 10° 
(dark-grey line), and A = 0.2 ± 0.1 K and a = 2.71 ± 0.03 
for regions with |b| > 20° (l ight-grey line). The noise level of 
the survey data, as given by lonas et al. (1998) , is also shown 
in Fig. H (dashed line). From the diagram it can be seen that the 
noise level is well below the signal, in the Z -range considered. 

The flattening of the spectrum at intermediate values of I 
is due to the presence of point sources in the map. In fact, the 
APS of a population of point sources randomly distributed in 
the sky is Poissonian, that is a spectrum with constant power 
at all I, As the fluctuation amplitude of the diffuse emission 
diminishes at small angular scales, the Poissonian component 
due to the point sources becomes dominant. 

Therefore, in order to analyse the angular power spectrum 
of the Galactic diffuse emission the point sources should be re- 
moved from the maps. Each point source in the map has a pro- 
file reflecting the point-spread function (PSF) of the observing 
beam. On a sufficiently resolved map one could attempt to re- 
move the sources by fitting the pixel temperatures with the PSF. 
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Fig. 1. Mollweide projection of the raw data from the Rhodes/HartRAO survey by Jonas et al. (1998), re-sampled into an 
HEALPix map. The temperature scale has been histogram-equalised. 




31,1 K 



Fig. 2. Histogram-equalised Mollweide projection of the Rhodes/HartRAO survey, after median filtering 



However at this low resolution several sources may be present the telescope beam size just over 3<7b. Median filtering, is an ef- 

in any given pixel, thus precluding this treatment. As a way of fective and easy-to-implement method to remove point sources 

suppressing the point sources from the data we have made use and suppress random instrumental noise (e.g. [Pratt 1991 ), with- 

of a median filter, that is, a convolution of the map with a me- out losing information at scales greater than three times the 
dian box of 9x9 pixels. A box of size 9 pixels corresponds to 
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Fig. 3. The angular power spectrum of the raw data from the 
Rhodes/HartRAO radio continuum survey at 2326 MHz for dif- 
ferent Galactic cut-offs. Three spectra are shown: all observed 
sky (black line), galactic regions with |6| > 10° (dark-grey 
line) and galactic regions with |6| > 20° (light-grey line). The 
survey noise level is also shown (dashed line). 



beam FWHM. The image of the Rhodes/HartRAO survey after 
being filtered is displayed in Fig. |[ 

The angular power spectra of the Rhodes/HartRAO survey 
for different latitude cut-off, after applying median filtering, are 
plotted in Fig. [|. For all the observed data, up to I — 100, the 
best-fit power-law to the APS has a amplitude A = 0.9 ± 0.2 K 
and spectral index a = 2.43 ±0.01. For CMB studies one 
is particularly interested in the properties of radio emission at 
high Galactic latitudes because this is where the Galactic sig- 
nal is weaker and where sensitive CMB measurements can be 
made. At this frequency the high Galactic latitude diffuse emis- 
sion is dominated by synchrotron radiation. Indeed the power 
spectrum of the emission appears to steepen if we mask out the 
Galactic plane. For a Galactic cut with |6 cu t| = 10°, a linear 
least-square fit to the APS up to I = 100 (dark-grey line) gives 
A = 0.3 ± 0.2 K and a — 2.81 ± 0.06; for \b cut \ = 20° (light- 
grey line) one obtains A = 0.3 ± 0.2 K and a = 2.92 ± 0.07. 
The quoted error takes into account an assessment of the un- 
certainties introduced by the partial sky coverage, the pres- 
ence of a Galactic cut and median filtering. We have evaluated 
this error by performing various simulations, as described in 
Sect. H Note from Fig. [| how the rapidly falling APS at high 
Galactic latitude becomes dominated by the aliasing noise from 
the equatorial cut off at I > 200: suppression of the sky signal 
by the beam-window function is no longer visible. 



Fig. 4. The angular power spectrum of the data from the 
Rhodes/HartRAO survey at 2326 MHz, after the point sources 
have been removed using a median filter. Line shading as in 
previous figure. 



2. 1. Point-source contribution 

The difference map of the raw data and the median-filtered 
data gives the map of the point sources plus the high-spatial- 
frequency instrumental noise. We have computed the APS of 
this difference map at high Galactic latitude and compared it 
with the APS from the sour ce counts expected at 2.3 GHz. 

Ibttolatti et al. (1998) report the point-source differential 
counts at 1.4, 5 and 8.44 GHz, for flux ranges between 10~ 5 



and 10 Jy. It can be seen from Fig.l of rpffolatti et al. (1998) 



that a differential source count distribution of the form 

dN/dS = 150S" 2 - 5 sr^Jy" 1 



(3) 



can be taken as a reasonable approximation of the source 
counts at 2.3 GHz in the flux range 0.1 — 10 Jy. Sources with 
lower flux levels will not contribute significantly to the fluctu- 
ation of antenna temperature in a 20 arcmin beam. The reason 
being that there are so many of these sources in a 20 arcmin 
beam that changing pointing does not cause significant changes 
to the average signal (Franceschini et al. 1989). In addition, 



the estimated receiver noise level for the Rhodes/HartRAO sur- 
vey is 25 mK and the measured Point Source Sensitivity (PSS) 
for the survey is 9.8 Jy/K. This means that point sources with 
fluxes lower than 0.2 Jy at 2.3 GHz are below the map noise 
level. 

We produced 10 realisations of 40,000 point sources with 
source counts of the form given in Eq. |[ in the flux range 
0.1 — 100 Jy, randomly distributed in the sky[| We convolved 

2 In order to assess whether 10 realisations were enough to charac- 
terize the population of point sources we generated another 20 reali- 
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the realisations with a Gaussian beam with 20 arcmin FWHM[] 
and computed the APS for each one. The APS of a randomly 
distributed population of point sources is given by 



C, 



dN 

Is 



S 2 dS 



(4) 



where S c is the flux cut ( Iegmark & Efstathiou 1996 ). 
Therefore the amplitude of the APS of different realisations of 
point-sources with the same differential source counts distribu- 
tion will vary slightly accordingly to the fluxes of the brightest 
sources present. 

In Fig. U we show the APS of the sources extracted from 
the Rhodes/HartRAO survey by median filtering for sky re- 
gions with \b\ > 20° (black line). The results can be compared 
with the mean APS of the 10 realisations of sources with dif- 
ferential counts distribution given by Eq. |] (continuous grey 
line). In the figure the standard deviation (SD) of the ampli- 
tude of the APS of the 10 realisations of point sources is also 
reported (dashed grey line). The APS of the point sources ex- 
tracted from the Rhodes/HartRAO survey is consistent with a 
population of point sources having differential source counts of 
the form given in randomly distributed in the sky. 

Comparing angular power spectra is not the ideal way of 
comparing source count distributions. As can be seen from 
Eq. [| the C; coefficients are given by the integral of the dis- 
tribution up to a flux-cutoff, therefore different source count 
distributions can give angular power spectra of the same am- 
plitude. However the APS of the point-source distribution is 
the quantity used to predict the level of contamination to CMB 
anisotropy c oming from extra-g alactic point sources (Toffolatti 
et al. 1998; |Hobson et al. 1999|). This result indicates that the 



basic assumptions on the spatial and flux distribution of extra- 
galactic radio sources used in CMB studies reproduce correctly 
the observations at 2.3 GHz made with a 20-arcmin FWHM 
beam. 

3. Simulations 

The incomplete sky coverage of the Rhodes/HartRAO data, 
point source suppression by median filtering and the presence 
of a Galactic cut affect the accuracy with which the angular 
power spectrum at high Galactic latitude may be recovered. In 
order to assess the impact of these factors on the estimates of 
the APS of diffuse radio emission we have performed a series 
of simulations and reproduced the processing procedure used 
on the real data. 

We have produced two sets of simulations: one for an input 
APS with a spectral index of a = 2 and the other for an input 
APS with a = 3. For each set we generated 10 different real- 
isations of a Gaussian random field on the sphere having that 
input APS. For this purpose we used the SYNFAST software 
(also part of the HEALPix software package). All the realisa- 
tions were generated at a resolution of 20 arcmin (FWHM) and 

sations and found that the mean of all the 30 realisations was within 
one standard deviation of the 3 subsamples of 10 realisations each 

3 The convolution was performed in the spatial-domain, to prevent 
any aliasing effect due to the point-like nature of the source signal. 
The convolving Gaussian beam was truncated at 5a 




■ Rhodes/HartRAO survey |b| > 20° 

■ Mean of simulated APS 

2-SD of the amplitude of simulated APS 



Fig. 5. APS of the point sources extracted from the 
Rhodes/HartRAO survey by subtracting the median filtered 
data from the raw data (black solid-line), away from the galac- 
tic plane. The mean of the APS of 10 realisations of a random 
population of point sources with differential source counts of 
the form given in Eq. 3 is also reported (grey solid-line). The 
grey dashed lines give 2 SD (standard deviation) of the ampli- 
tude of the APS from the simulations. Model and observations 
are consistent with each other. 



sampled into HEALPix maps with nside = 5 12, as in the case of 
the Rhodes/HartRAO map. From these we computed the APS 
for three different conditions: total sky coverage, Galactic plane 
cutoff with |& cu t| = 10° and partial sky coverage plus Galactic 
cutoff ( | £»cut | = 10°). For an input APS with a = 3 we also 
considered the case in which a population of spatially random 
point sources is added onto the simulated diffuse emission. 



3. 1. Diffuse emission only 

In Fig. ||the means of the derived APS from 10 different re- 
alisations of a spectrum with spectral index a = 2 are shown. 
The curves are for the three different cases considered, that is 
deriving the APS from full sky maps of the realisations, from 
maps with a Galactic cutoff with |6 cu t| = 10° and from maps 
with coverage identical to that of the Rhodes/HartRAO survey 
plus a Galactic cutoff. With a Galactic cutoff both the usual 
definition for the APS and the Modified APS definition have 
been used. As seen from Fig. [g, in this case, the incomplete 
sky coverage and the presence of a Galactic cutoff do not affect 
the recovery of the input APS in a substantial way, and there is 
no appreciable difference between the standard APS definition 
and the Modified definition. 
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Fig. 6. Mean recovery of the APS from 10 different realisations of an input APS with spectral index 2 (dashed light-grey line), for 
the three different cases considered: full sky coverage, galactic cut at \b\ = 10° and Rhodes/HartRAO sky coverage plus galactic 
cut at \b\ — 10°. The dark-grey lines refer to the standard definition for the APS, the black lines to the Modified APS. 
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Fig. 7. The mean of the spectra recovered from 10 different realisations of an input APS with spectral index —3. The curves are 
as in Fig. 6. 



In order to quantify the accuracy of the recovery we have 
computed the best-fit power law of the spectra obtained from 
the 10 realisations, in the three different cases of sky coverage. 
The amplitude and spectral index of these best-fit power laws 
can be compared with the true amplitude and spectral index 
of the APS used as input for the 10 realisations. The standard 
deviation of the derived amplitudes and spectral indices with 
respect to the true (input) values provides an estimate of the 
standard error on the recovered amplitude and spectral index. 
The results are summarised in Table [j]. For a full sky cover- 
age the uncertainty in the measured spectral index over the l- 
range 2-200 is of about 2% as can be seen simply by taking 
the standard deviation over the mean. An incomplete sky cov- 
erage such as the one of the Rhodes/HartRAO survey plus the 
introduction of a Galactic latitude cutoff brings this uncertainty 
up to 4%. The error on the amplitude of the power law is not 
significantly affected by a diminished sky coverage. 

The same cases of sky coverage have been considered also 
for an input APS with spectral index a = 3. Fig. [7] shows the 
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Table 1. Mean and standard deviation of the recovered ampli- 
tude and spectral index from 10 different sky realisations of 
an input APS with amplitude A = 1.0 K and spectral index 
a = 2 (top panel) and a = 3 (lower panel). Three different 
cases have been considered: full sky coverage, Galactic cut at 
|&cut| = 10° and Rhodes/HartRAO sky coverage plus Galactic 
cut. I G [2,200]. 



mean of the recovered spectra from 10 different sky realisa- 
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tions of the input APS. In this case the derived spectra are more 
strongly affected by a Galactic cutoff and an incomplete sky 
coverage. The modulation of the APS caused by introducing 
the Galactic plane cutoff is clearly visible in the figure (middle 
panel). The figure also shows how the use of the Modified APS 
definition helps in reducing such aliasing. 

The convolution of the field Spherical Transform by the 
Spherical Transform of the sky coverage window (which is 
undersampled) introduces a coupling between the a; TO coef- 
ficients of the field, which is stronger for fields with steeper 
APS, simply because the gradient in I is stronger. This is the 
reason why an input APS with a = 2 is barely affected by the 
incomplete sky coverage, while this has such a strong impact 
on the measure of an input APS with a = 3. In practice the 
coupling introduced by the sky-coverage window means that 
no APS with spectral index steeper than —3.5 can be recovered 
in the presence an equatorial cutoff or large areas of missing 
data. 

However, despite the aliasing noise introduced by the 
Galactic plane cutoff and the incomplete sky coverage, one can 
still recover the spectral index of an input APS with a ~ 3. The 
parameters resulting from least-square fits to the Modified APS 
obtained from the 10 realisations are reported also in Table [j]. 
The percentage error on the recovery of the input spectral index 
from a full-sky observation, is in this case about 1.5% and does 
not increase substantially with decreasing sky coverage. 



(diffuse emission plus sources) are respectively a = 2.49 ± 
0.02 and a = 2.79 ± 0.02^, and after median filtering one 
obtains a = 3.04 ± 0.07 and a = 2.96 ± 0.07, which are good 
recoveries of the true spectral index of the simulated diffuse 
emission. 
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3.2. Diffuse emission and point sources 



In Sect 2. 1 we have described how we generated 10 realisations 
of a spatially random population of point sources with differ- 
ential source counts of the form given in Eq. ^| From Toffolatti 
et al. (1998) one can see how this expression of the differential 
source counts approximate quite well the extra-galactic source 
counts that are expected at 2.3 GHz, in the flux range 0.1 — 10 
Jy. Sources with flux lower than 0. 1 Jy do not contribute signif- 
icantly to the temperature fluctuations in the Rhodes/HartRAO 
survey. The 10 different realisations of point sources were then 
added onto 10 different realisations of an input APS with a 
spectral index of a = 3, normalised to the amplitude of the 
APS of the Rhodes/HartRAO survey measured for a Galactic 
cutoff with |b cut | = 10°. Each resulting map, which mimics 
the superposition of diffuse emission with a Poissonian pop- 
ulation of point sources was then convolved with a median 
filter of size 9x9. The APS and the Modified APS of the fil- 
tered map was then computed. This APS represents a recovery 
of the spectrum of the diffuse component and can be directly 
compared with the input spectrum of the diffuse component to 
assess the systematic error introduced by median filtering the 
map. 

In Fig. U the angular power spectra of two simulated skies 
are shown on the same scale as the data in Fig. |3|. As in the 
case of the real data the flattening of the spectrum at high I due 
to the Poissonian APS of the point source component is clearly 
recognizable. These 2 particular cases of simulated skies mimic 
quite closely the spectra of the raw data in Fig. [}J In these two 
cases the fitted-spectral index up to I = 100 of the input maps 



Fig. 8. APS of 2 simulated skies containing a diffuse compo- 
nent with power law spectrum with a = 3 and a population 
of point source of the form given in Eq. 3. Note the similarity 
with the APS of the raw data at high Galactic latitude shown 
in Fig. 3. In these two particular cases the spectral index up to 
I = 100 are a = 2.49 ±0.02 and a = 2.79 ±0.02 respectively, 
after median filtering the maps one obtains a — 3.04 ± 0.07 
and a = 2.96 ±0.07. 



Fig. ^| shows the Mean APS of the diffuse component re- 
covered by median filtering in the three different cases of sky 
coverage considered. The median filtering suppression of the 
high frequency signal starts to affect the APS at I ~ 100. Up 
to this multipole order the filtering process allows the correct 
spectral index of the input APS to be recovered, with an error 
lower than 3%, as taken from Table |[ The average spectral in- 
dex of the spectra of the 10 maps before median filtering is 2.8, 
with a standard deviation of 0.1 for I 6 [2, 1001. 



4. Discussion 

The analysis of the Rhodes/HartRAO data shows that the global 
APS of synchrotron emission at high Galactic latitudes (|6| > 
20°) is well approximated by a power law spectrum having 
spectral index a = 2.92 ± 0.07. This result at 2.3 GHz is con- 



sistent with the results of Iegmark & Efstathiou (1996) and 



the fit error is given in this case 
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Fig. 9. Mean of the spectra from 10 different realisations of an input APS with spectral index a = 3 (dashed light-grey line) plus 
a population of point sources with differential source counts of the form given in Eq. 3. The maps were first processed by median 
filter to remove the point source components. The three different cases of sky coverage considered are shown. The dark-grey 
lines refer to the standard definition for the APS, the black lines to the Modified APS. 
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Table 2. Mean and standard deviation of the 
tude and spectral index of the APS of the diffu 
from 10 different sky realisations. Each realisa ion 
of a realisation of an input APS with amplitude 
and spectral index a = 3 plus a realisation of a 
dom population of point sources with source 
given in Eq. 3. The mean spectral index of the A 
sky unprocessed maps was 2.8 ± 0.1. / G [2, 100 



! coui ts 

l 5 



foi nd 



' pate h 



Bouchet et al. (1996) who found a spectral inde) 
up to I ~ 100 from the analysis of square patch 
MHz survey. Bouchet & Gispert (1999) also 
I > 100 from the analysis of patches of the 1420 

Lasenby (I99'7) found a ~ 2 by analysing 
and the 1420-MHz surveys in the "Tenerife 
is a 60 x 10 degree patch of the northern hemi: 
Galactic latitude. This could be a local property 
of the sky. Alternatively, it could be due to the 
sources: Lasenby did not attempt to subtract poin : 
vate communication) and this may have the eff 
ing the spectrum. From Fig. [| one can see that 
source subtraction the APS of the synchrotron 
slope of approximately —2 in the /-range 50-15C 
/-range that can be probed with a patch of sky of 
"Tenerife patch", at a resolution of ~ 1°. 

A spectral index of a ~ 3 has also been d 



APS of Galactic dust emission ( 3autier et al. 1 )92 



& Efstathiou 1996; |Wright 1998| ) and this is 
measured f or the APS of single velo city neutr 
the Galaxy ( Crovisier & Dickey 1983 ). However 
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(1998) report a flatter slope of a = 2.5 for dust emission from 
the co nbined DIRBE and IRAS dust maps, at high Galactic 
latitudes. 

simulations we have performed showed that incom- 
coverage and a Galactic cut do not introduce signif- 
ors when deriving the spectrum of an input field with 
~ 2 . This is reassuring since a ~ 2 is approximately the 
1 index of the APS of the CMB. On the other hand the 
sky coverage of the Rhodes/HartRAO survey and a 
c cut seriously degrade the quality of the spectrum that 
derived from an input field with a spectral index steeper 
Therefore, while spherical harmonic decomposition is a 
i tool to analyse the statistical properties of CMB emis- 
incomplete sky coverage, it is not ideal when one 
analyse fields whose spectra have index a > 3. This 
hts the need to refine the analysis tools in view of the 
vations of the Galactic diffuse emission that will be 
out by Planck and MAP satellites. 

one limit of our simulations is that the synthesized 
emission is a realisation of a random Gaussian field 
given APS while the Galactic signal is not a spatially 
field. The spectral index measured from the map of 
c radio emission changes significantly as the brightest 
of the Galactic plane are gradually masked off. This 
happen when one applies the same type of mask- 
a random field. We are planning to study the level of 
of the Galactic signal from the Rhodes/HartRAO 
data in a future work. 

it is apparent from Fig. |2|, the Rhodes/HartRAO survey 
affetted by a certain level of striping, which reflects base- 
stabilities due to ground spillover and weather effects 



(Jonas et al. 1998). This is more visible around the South Pole 
region where the scans necessarily traversed low elevation an- 
gles. This striping, with a typical angular scale of ~ 5°, will 
introdi ce some additional power in the spectrum at multipole 
order , ~ 40. However, as apparent from the map, even at 
high-C alactic latitudes, the amplitude of these striations are on 
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the average lower than the observed variation of sky brightness 
temperature, therefore we expect striations to be only a minor 
contribution to the overall power measured around Z ~ 40. The 
fact that no significant change of slope is observed in the APS 
after point-source suppression for I < 100 (Fig. ||), is another 
indication that striping makes a minor contribution to the total 
power at a given angular scale. If striping does make a ma- 
jor contribution, then its effect is being disguised by a steeply 
falling power from the Galactic diffuse emission. So, for the 
purpose of estimating the contribution of synchrotron emission 
relatively to the CMB anisotropy, a spectral index of a = 2.9 
for the synchrotron APS is a conservative estimate. 

One can use the APS derived at 2.3 GHz to infer the level of 
contamination of the CMB signal at the different angular scales 
at 30 GHz: this is the lowest frequency channel of Planck and 
the one for which the highest level of synchrotron con tamina- 
tions are expected. On the 7° scale, Kogut et al. (1996) derived 
an upper limit of 1 1 [iK to the temperature fluctuations due to 
synchrotron emission in the COBE DMR channel at 3 1 .5 GHz 
for Galactic latitudes with |6| > 20°. If we combine this limit 
with our measure of a = 2.92 ± 0.07 for the spectral index of 
the APS we derive an upper limit of A30GHZ = 120 /iK for the 
amplitude of synchrotron emission . This upper-limit to syn- 
chrotron emission is shown against a CMB spectrum in Fig. |l(| 
as a function of the multipole order /. The APS for the CMB has 
been computed assuming a flat inflationary model with purely 
scalar scale-invariant fluctuations, vacuum energy S1a = 0.43, 
cold dark matter density /i 2 f2 c d m = 0.20 and baryon density 
h 2 Q,\y = 0.03, (where h = 0.63 is the Hubble parameter), as 
derived by [Tegmark & Zaldarriaga (2000) from the latest CMB 
measurements ( ge Bernardis et al. 2000| ; |Winant et al. 2000| ). 
From Fig. |l(| it is apparent that synchrotron emission can be 
a significant contribution to the sky anisotropy at 30 GHz at 
scales larger than ~ 20°, but it is an order of magnitude weaker 
than the cosmological signal at degree angular scales where 
the first "acoustic" peak of the CMB is observed. Indeed from 
COBE observations we know that for the quadrupole (1=2), 
Galactic emi ssion is comparable in amplitude to the anisotropy 
of the CMB ( fcennett et al. 1996j ). 

If — (3 is the spectral index of the frequency dependency 
of A, then -A30GHZ = 120 /iK combined with the value 
^2.3GHz = 0.3 ± 0.2 K measured at 2.3 GHz provides a lower 
limit on the frequency spectral index of f3 > 2.9, between 2.3 
and 31.5 GHz. At large angular scales (18°), Platania et al. 
(1998) derive an average value for the spectral index of the 
synchrotron frequency dependency between 1 and 7.5 GHz of 
j3 = 2.81 ± 0.16. The spectral index is expected to steepen 
at frequencies higher than 10 GHz, since the observed en- 
ergy spectrum of cosmic ray electrons steepens at energy of 
- 10 - 20 GeV ( |Banday & Wolfendale 1990| ; Banday et al. 
1991). 

In order to bracket the synchrotron contribution to the sky 
fluctuations at 30 GHz one can combine the lowest (ler) bound- 
ary of our measure for the amplitude, ^2.3GHz = 0.1 K, 
with the lowest boundary of the frequency spectral index de- 
rived by platania et al. (1998) , that is f3 = 3: one obtains 
^30GHz = 40 /iK. This lower limit on the synchrotron con- 



tribution is also shown in Fig. [T(] as a function of the multipole 
order. 

Tegmark et al. (2000) give three possible estimates of A 
for Galactic synchrotron emission at 19 GHz and label them as 
optimistic, pessimistic and "middle-of-road". Their pessimistic 
value of 192 /iK at 19 GHz with a frequency spectral index of 
(3 = 2.6, corresponds to a ^oghz ~ 60 /iK which is within 
the range of our estimate, but somewhat on the lower bound. 
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1 000 



Fig. 10. The estimated range for the APS of synchrotron emis- 
sion at 30 GHz (grey line) can be compared here with the APS 
of the CMB (black line), as derived with the latest estimate for 
the cosmological parameters (see text). In the Z-range where 
it has been extrapolated the APS of synchrotron emission is 
shown as a dashed line. The square point is the upper limit 
to synchrotron temperature fluctuations derived from COBE 
DMR at 31.5, for the 7° angular scale. 



5. Summary and Conclusions 

In order to analyse the global angular power spec- 
trum of Galactic radio emission we have re-sampled the 
Rhodes/HartRAO Galactic survey at 2326 MHz into a 
HEALPix map with nside= 512. The point sources were re- 
moved from the map by median filtering. We find that a best-fit 
to the angular power spectrum of the entire Rhodes/HartRAO 
survey with a power-law spectrum Ci oc l~ a , gives a = 
2.43 ± 0.05 in the I range 1 - 100. At high Galactic latitudes, 
|6| > 20°, where CMB observations are carried out, the APS 
of the radio emission steepens and it is well described by a 
power-law with spectral index a = 2.92 ± 0.07 and amplitude 
^2.3GHz = 0.3 ± 0.2 K, up to I = 100. At these latitudes radio 
emission is dominated by synchrotron emission. By combining 
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our results with the upper limit of 1 1 /iK on the synchrotron 
contribution to temperature fluctuations at 3 1 .5 GHz, obtained 
by COBE-DMR for the 7-degree angular scale, we derive an 
upper limit of ^4.3oghz = 120 /iK to the amplitude of syn- 
chrotron emission at high Galactic latidue. This in turn implies 
a lower limit of (3 > 2.9 on the spectral index of the frequency 
dependency of A. 

With an APS index of a = 2.9 and an amplitude A30GHZ < 
120 /iK, at 30 GHz synchrotron emission could be a contami- 
nant of CMB observations at large angular scales (/ < 10), but 
will not obstruct a clean view of the CMB anisotropy at degree 
angular scales. 

A map of point sources from the Rhodes/HartRAO survey 
was also obtained. The APS at high Galactic latitude was com- 
pared with the prediction derived from extra-galactic source 
counts, as interpolated between 1.4 and 5 GHz, under the as- 
sumption of a purely random sky distribution. The same source 
counts and the randomness assumption are the bases of the 
studies that predict the level of temperature fluctuation to be 
expec ted from extra-galactic point sources in CMB obs erva- 
tions dFranceschini et al. (1989| [Toffolatti et al. (1998)| ). We 
conclude that at 2.3 GHz and for a resolution of 20 arcmin, 
predictions of the point source contribution and observations 
are consistent with each other. 

The simulations that we have performed to evaluate the 
accuracy of the results show that the spectra derived from 
the Rhodes/HartRAO survey in the /-range 2-100 are reliable. 
Incomplete sky coverage does not compromise the accuracy 
with which one can derive the APS of the field by spherical har- 
monic decompoistion, if this field has a spectrum with spectral 
index a < 3. However we have also found that the quality of 
the spectrum derived degrades significantly for spectral index 
a ^ 3, if large areas of the sky are missing. 

The simulations showed that median filtering can be used 
to effectively remove point sources and derive the APS of the 
diffuse component up to I = 100, for a beam FWHM of 20 
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